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Unusual dependence of vortex core states on the superconducting gap 
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We present a scanning tunneling spectroscopy study on quasiparticle states in vortex cores in 
Bi2Sr2CaCu208+i . The energy of the observed vortex core states shows an approximately linear 
scaling with the superconducting gap in the region just outside the core. This clearly distinguishes 
them from conventional localized core states, and is a signature of the mechanism responsible for 
their discrete appearance in high-temperature superconductors. The energy scaling of the vortex 
core states also suggests a common nature of vortex cores in Bi2Sr2CaCu20g+i and YBa2Cu307_i. 
Finally, the observed vortex core states do not show any dependence on the applied magnetic field 
in the range from 1 to 6 T. 

PACS numbers: 74.50.+r, 74.60.Ec, 74.72.Hs 



In conventional, s-wave superconductors, the suppres- 
sion of the order parameter in a vortex core creates a po- 
tential well for low-energy quasiparticles, leading to the 
formation of localized states On the contrary, if 

the superconducting order parameter has nodes in it - as 
for dx2-y2 symmetry in high-temperature superconduc- 
tors (HTS) - one expects the low-energy quasiparticles 
states in a vortex to be extended along the nodes of the 
gap function, so to be delocalized. This would result in 
a broad peak at the Fermi level in the quasiparticle lo- 
cal density of states (LDOS) of a vortex core In 
this light, the observation of discrete vortex core states 
in YBaaCusOy-s (YBCO) § and BiaSraCaCusOg+s 
(BSCCO) [|||] by scanning tunneling spectroscopy (STS) 
has come as a complete surprise. As a result, the nature 
of these states has been subject of increasing theoretical 
study leading to a range of possible scenarios 

for explaining the experimental data. The need for un- 
derstanding of the electronic structure of the vortices in 
HTSs has become even more pressing because of the an- 
tiferromagnetic fluctuations recently observed in vortex 
cores in La2„a;Sr3,Cu04 |17|. 



To our knowledge, the most direct way to access the 
electronic structure of a vortex core is by using a scanning 
tunneling microscope (STM). In a typical experimental 
set-up for studying vortex cores in HTSs (as well as in this 
study), the STM tip and tunneling direction are perpen- 
dicular to the CuOa planes. The quasiparticle excitation 
spectrum then follows from the dl/dV tunneling spectra 
p8[ . Experiments on HTSs have revealed the following 
characteristics of vortex cores. In YBCO the vortex core 
states appear as two clearly distinct, sub-gap excitations, 
which do not disperse on moving out of the vortex core, 
but rather transform into weak shoulders in the super- 
conducting spectra (and sometimes also observed at zero 



magnetic field) 0. In BSCCO the vortex core spectra 
reveal a remarkable resemblance to the pseudogap spec- 
tra observed above the critical temperature Tc [^. In 
addition, vortex core states appear as weak shoulders in 
this pseudogap, do not change in energy as a function 
of increasing distance from the vortex core |^ (contrary 
to the vortex core states in conventional superconduc- 
tors Q), and decay over a characteristic length scale of 
about 20 A |9|. Though quite irregular shapes can be 
observed due to vortex motion , the cores do no show 
any sign of the four-fold symmetry that may be expected 
for d-wave superconductors. 

In this letter, we present new data on vortex cores 
in BSCCO, obtained with a low-temperature STM. Our 
main results can be summarized as follows. First, the 
energy of the vortex core states scales with the super- 
conducting gap outside the core, clearly distinguishing 
these states from localized vortex states {E oc A'^/Ep) 
in conventional superconductors . Since in BSCCO the 
superconducting gap scales with the oxygen doping level 
pl|-p6[ , this directly gives the doping dependence of the 
vortex states as well. Second, the vortex core spectra 
do not show any significant dependence on the external 
magnetic field over a range from 1 to 6 T, which ques- 
tions field-dependent scenarios for explaining the vortex 
core states. 

A systematic study of vortex core spectra requires, 
apart from sufficient instrumental resolution, relatively 
high sample homogeneity. More precisely, in order to 
compare spectra inside a vortex core to those in the 
nearby superconducting region, (zero-field) spectral re- 
producibility over at least 100 A is necessary. This con- 
dition can be met in BSCCO single crystals with tran- 
sition widths ATc < 1 K (as measured by AC suscepti- 
bility). Most of the data presented here were obtained. 
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with magnetic field and tunneling direction perpendicu- 
lar to the Cu02 planes, on two different overdoped sam- 
ples: (i) A sample with 77.7 K (AT^ = 0.4 K), 
cooled down at 6 T, and after measurements at 6 T the 
field was reduced (at low temperature) to 4 and 1 T. 
(ii) A sample with = 77 K (AT^ = 1 K), zero-field 
cooled, measured at fields of subsequently 0, 6, 2, and 
T. For the latter sample, lines of spectra in and near 
vortex cores, as well as vortex images as a function of field 
can be found in Refs. |P,pO|. Some additional data were 
taken on two other overdoped samples with Tc = 76.1 K 
(ATc = 0.3 K), and on an optimally doped sample with 
Tc = 87.4 K (ATc = 1.0 K). Samples were cleaved in 
ultra-high vacuum environment (10~^ mbar), shortly be- 
fore cooling down STM and sample to 4.2 K in exchange 
gas 10~^ mbar helium). All measurements in this 
study were taken at 4.2 K. For further experimental de- 
tails we refer to Refs. |p7|-p9|. 

The quasiparticle LDOS at an energy E — eV on the 
surface of the sample was obtained by measuring the dif- 
ferential tunneling conductance d//dy at sample bias V 
as a function of position, using a lock-in technique. In ho- 
mogeneous BSCCO samples, clearest vortex images were 
found by measuring the conductance at V = — Ap/e, 
where Ap corresponded to the energy of the supercon- 
ducting coherence peaks. We thus obtained maps of spots 
where superconductivity was suppressed. A comparison 
of the number of these spots per area to the magnetic fiux 
for different fields, as well as their disappearance at zero 
field, justified the identification as vortex cores |19,M. 




X [nm] 100 

FIG. 1. Image of vortex cores at 6 T (field cooled mea- 
surement), taken by measuring the conductivity dl/dV at 
-30 mV, and normalized on its value at zero bias. 

In Fig. |l| we show a map of vortex cores at 6 T. As can 
be verified directly, the number of vortices (30 ± 2) in 
Fig. corresponds to the average flux crossing this area 
(29<I>o) at 6 T. In general, the density of vortices scales 
with the magnetic field as one should expect. We do not 
observe any systematic change in the size and shape of 
the vortex cores for the different applied fields. 

dl/dV spectra were taken in and around several of 
the thus imaged vortex cores. Spectra along a 200 A 
trace through a vortex core are shown in Fig. ^a). The 
spectra discussed hereafter were obtained by averaging 



the spectra just outside the core [above and below in 
Fig. ^(a)] and the spectra well inside the core [in the 
middle of Fig. |^(a)], for the superconducting spectra and 
the typical vortex core spectra respectively. These av- 
eraged spectra can be found in Fig. ^(b-d), for a vortex 
core at 1 T and for one at 6 T in BSCCO, as well as 
for a vortex core in YBCO 0|). Compared to the su- 
perconducting state, the LDOS in the vortex cores seems 
considerably reduced, with loss of spectral weight near 
the Fermi level. The characteristics of the vortex spectra 
in BSCCO, asymmetric pseudogap with weak shoulders 
at low bias, do not show any dependence on the field 
strength. We neither observe, within the experimental 
resolution, any change (due to the magnetic field) in the 
spectra outside the vortex core. The spectra taken at dif- 
ferent field strengths are in fact remarkably similar, both 
those inside, and those around vortices. 
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FIG. 2. (a) Spectra along a 200 A trace through a vor- 
tex core in BSCCO at 1 T. The vortex core is about halfway 
the trace. The spectra have been ofi'set for clarity, (b) Aver- 
aged spectrum in the centre of this core, and in its immediate 
vicinity (dashed line). The latter spectrum may appear a bit 
broad, since it is the average of several spectra with peaks at 
slightly difi'erent energies, (c) For a vortex core in BSCCO, at 
6 T. (d) For a vortex core in YBCO at 6 T. Note the diflterent 
scales for the bias voltage. 



In order to study the dependence of the vortex core 
spectra on the superconducting energy gap (taken as the 
energy of the coherence peaks Ap), we measured vor- 
tex cores on different samples, and also made use of a 
certain inhomogeneity (on a scale larger than 100 A) of 
some of the samples which we would characterize as mod- 
erately homogeneous. (In fact large-scale homogeneity 
has also been obtained, using an optimum oxygenation 
procedure.) Ap has been determined from the average 
of spectra around the vortex core, and its precision de- 
pends on the inhomogeneity of the superconducting order 
parameter in immediate neighborhood of the core. 
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For a precise and consistent determination of the ener- 
gies of the core states, the vortex core spectra were fitted 
with a fifth order polynomial over a range 0.6Ap < eV < 
1.5Ap, with the fits forced to go through the zero-bias 
conductance. These fits of the pseudogap background 
were subtracted from the vortex spectra, thus leaving 
the excess spectral weight related to vortex core states. 
For the vast majority of vortices, the core state ener- 
gies as determined from the fitted and subtracted data 
fall within the error bars of the core state energies de- 
termined directly from the raw data (see the arrows in 
Fig. |). The latter are less precise because the core states 
are partly hidden in the slope of the pseudogap |^ . The 
results have been checked for robustness against varia- 
tion of the fitting parameters and range, and have been 
plotted on a normalized energy scale in Fig. H. 




FIG. 3. Spectra of vortex core states obtained by subtract- 
ing a fitted background (see text) from the vortex core spec- 
tra, plotted with an energy scale normalized on the energy 
of the coherence peak Ap of spectra near the core. There is 
no significant variation in the spectra for different magnetic 
fields and A„. 



These data confirm the independence of the vortex 
states on the magnetic field, and also show that their en- 
ergy Ecore directly scales with the superconducting gap. 
This scaling comes out even more clearly when Ecore is 
plotted as a function of Ap, for all measured vortices 
in different samples and at different fields, see Fig. ^. 
The dependence on Ap is certainly not quadratic, as one 
would expect for conventional localized states in vortex 
cores [||. It is roughly linear, passing through the ori- 
gin, and with a slope of 0.28. We have checked that the 
energies extracted from the raw data, with a bit more 
scattering of the data and a slope of approximately 0.30, 
do show the same behavior. This linear behavior is an im- 
portant result, since it proves that the vortex core states 
in BSCCO do not correspond to conventional localized 
states. 

Since doping dependent tunneling experiments on 



BSCCO show a roughly linear dependence of Ap on oxy- 
gen doping pl|-p6t, we can directly read the horizontal 
scale as the (local) hole doping of the BSCCO samples, 
going from overdoped to underdoped. This indicates that 
the electronic nature of the vortex cores remains the same 
for a considerable doping range. 
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FIG. 4. The energy of vortex core states Ecore as a func- 
tion of the superconducting gap (the energy of the coherence 
peaks) Ap. Plotted are data for different magnetic fields, and 
different BSCCO single crystals. At Ap = 17 meV the value 
extracted from YBCO data. The straight line is a linear fit to 
all data points, forced to go through the origin, with a slope 
of 0.28. An unrestricted linear fit would cut the vertical axis 
at about 1 meV, the latter value also being the uncertainty 
of the offset. 

In Fig. I we have also included data obtained on 
YBCO, for comparison Q] . Though the vortex core states 
in YBCO appear much more pronounced than those in 
BSCCO (see Fig. their energy scale shows a remark- 
able similarity, strongly suggesting a common origin. It is 
tempting, but rather speculative, to compare this to the 
energy scale of the antiferromagnetic fluctuations with 
an energy of 3 4 meV, for a spin gap of 6.7 meV, as 
observed in La2_i;Sr^Cu04 [0. 

Our results, especially the linear scaling of Ecore with 
Ap, rule out conventional localized states as an explana- 
tion for the vortex core states observed in HTSs (they 
would neither be consistent with the lack of dispersion 
as a function of position, see Fig. Furthermore, they 
provide sufficient details to discuss other scenarios. 

First of all, in a pure d-wave superconductor the vor- 
tex core states would appear as a broad zero-bias con- 
ductance peak iQ-^, in clear disagreement with experi- 
ment. One may avoid this difficulty by supposing that 
the quasiparticle states from different vortex cores form 
bands, which would lead to the splitting of the zero-bias 
conductance peak into two different peaks. In that case, 
however, the energy gap between these peaks should de- 
pend on the magnetic field, which is inconsistent with 
the results presented above. 

A possible explanation for the observed vortex core 
states is the symmetry breaking of the order parameter 
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in the vicinity of a vortex core ||10|,|11[. This will lead 
to secondary (possibly imaginary) dxy or s components 
of the order parameter, effectively blocking the nodes of 
the dx'^-y^ gap function, and allowing localized states 
Though in a simple BCS d-wave superconductor these 
components are too small to influence the spectra 
they appear more important when the Coulomb interac- 
tion on the Cu sites is taken into account. Numerical 
calculations using the t — J model estimate their size - 
strongly doping dependent ~ to be between 5 and 30% of 
the Ax2_y2 order parameter ^3 15 1. Our results imply 
that such a secondary component would scale with the 
A2;2_j,2 gap over a wide doping range and, contrary to 
what one should expect |H,|l^, would be independent of 
the magnetic field between 1 and 6 T. 

In this context, one may ask whether the observed core 
states are directly related to the vortices, or are merely 
enhanced inside the cores. In the secondary order param- 
eter scheme, the latter possibility would mean that this 
secondary component is also present at zero field. In fact, 
low-energy features similar to vortex core states (though 
considerably less pronounced) have been observed in the 
gap of YBCO at zero field Q , and interpreted as the sig- 
nature of an additional s or dxy component in the order 
parameter in overdoped YBCO pO| . 

The energy scaling and field independence of the core 
states naturally suggest an other scenario, involving the 
pseudogap state. Keeping in mind that the pseudogap 
above Tc (and in the vortex cores) scales with the super- 
conducting gap as well 1 22 2^j2^,^ , one can speculate 
that the presence of the pseudogap leads to a splitting of 
the BCS d-wave zero bias conductance peak. This idea 
has been elaborated in two very recent calculations, us- 
ing the two-body Cooperon operator for modeling phase 
fluctuations and invoking spin-density wave order in 
the vortex cores Q. The pseudogap being much more 
dominant in BSCCO than in YBCO, it will suppress the 
low-energy core states in the former material much more 
than in the latter [Q, in agreement with our observa- 
tions (Fig. |). 

In conclusion, we have measured vortex core quasipar- 
ticle states as a function of the superconducting gap and 
the magnetic field. Their energy dependence on the su- 
perconducting gap near the vortex core is approximately 
linear, with a slope of about 0.28. The vortex core states 
do not show any dependence on the applied magnetic 
field in the range from 1 to 6 T. Our results suggest 
a common origin for vortex core states in YBCO and 
BSCCO. The linear energy scaling of these states with 
the superconducting gap, the lack of any dispersion as a 
function of distance from the vortex core centre, and the 
suppression of LDOS in the core compared to the super- 
conducting state further underline the non-BCS behavior 
of HTSs, and of their vortex cores in particular. 
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